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Nesting in a semimetal can lead to an excitonic insulator state with spontaneous coherence between 
conduction and valence bands and a gap for charged excitations. In this paper we present a theory of 
the ferromagnetic state that occurs when the density of electrons in the conduction band and holes 
in the valence band differ. We find an unexpectedly rich doping-field phase diagram and an unusual 
collective excitation spectrum that includes two gapless collective modes. We predict regions of 
doping and external field in which phase-separated condensates of electrons and holes with parallel 
spins and opposing spins coexist. 
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It has long been recognized that a semimetal or 
small-gap semiconductor could undergo a phase transi- 
tion to a state with spontaneous coherence between con- 
duction and valence bands and a gap for charged excita- 
tions. When described in terms of electrons in the con- 
duction band and holes in the valence band, this state 
can be regarded as an electron-hole pair (exciton) con- 
densate, and its mean-field theory description is very sim- 
ilar to the BCS mean-field theory of a superconductor. 
Excitonic insulator states are expected to occur only if 
conduction and valence bands are adequately nested. 

The excitonic insulator state is evidently rare, most 
likely J3j because accurate nesting is uncommon and be- 
cause electron and hole densities are usually not equal. 
The most compelling example of excitonic-insulator be- 
havior to date is that demonstrated by bilayer quantum 
Hall systems Q], where a charge gap develops in a state 
with spontaneous coherence || between electrons in dif- 
ferent layers. In this paper we consider doped excitonic 
insulators, i.e. systems in which the conduction-band 
electron density and valence-band hole density differ. 
Doping is unfavorable for electron-hole pairing, causing 
the excitonic insulator to evolve first into a metallic ex- 
citonic ferromagnet ||, and finally into a paramagnetic 
normal metal state. In this paper we address properties 
of the intermediate ferromagnetic state. 

This work is motivated in part by the recent discovery 
of weak ferromagnetism in lightly-doped divalent hexa- 
borides (?]]. These materials are ferromagnetic in spite 
of the absence of d- or /- partially filled orbitals, and 
have very high Curie temperatures T c > 600-K" that de- 
pend strongly on carrier concentration. Early electronic 
structure calculations had predicted that these materi- 
als would be semimetals or small gap semiconductors [|| , 
leading Zhitomirsky et al || to propose excitonic pairing 
as the source of their ferromagnetism, renewing inter- 
est in this unusual magnetic state jH||l|. Uncertainty 
has been generated, however, by recent experiments |r"^ ] 
which appear to imply that these materials are semicon- 
ducting with a ~ leV gap, in agreement with more re- 



cent GW approximation electronic structure calculations 
ifLsf , casting doubt on the excitonic insulator picture. In 
this Letter, we present an analysis of the collective ex- 
citation spectrum, and doping-chemical potential phase 
diagram of excitonic insulator ferromagnets, establish- 
ing features which could be used to convincingly identify 
these states on purely experimental grounds. In particu- 
lar, we demonstrate that doped excitonic condensates are 
unusual ferromagnets with a characteristic multi-branch 
collective mode structure. 

We start by considering a model Hamiltonian H = 
Ho + Hi where Ho includes both conduction and valence 
band energies and Hi includes all interactions. 

Ho = ^ [(eka - - hr CT:0 ./] 4. acr c ka(T / (1) 
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Here k denotes the quasimomentum, a the spin, a = c.v 
is a band label, fj, is the chemical potential and h = 
ig/isH here H is an external magnetic field. We start by 
including only long-range Coulomb interactions, which 
leads to a model with interaction vertices that conserve 
band labels and spins: 

Ul = ^n H^( q ) C k+qa ff C L-qaV' C k'aV'C ka(T (2) 

kk'qerer' aa' 

where f2 is the crystal volume. This interaction model is 
known as the dominant term approximation and is 
accurate for low-carrier densities and weak Coulomb in- 
teractions. In this model charge and spin are conserved 
separately in each band, leading to a SU (2) x SU(2) sym- 
metry [ fl0| that will figure prominently in our discussion. 

For simplicity we consider the case of delta-function 
interactions, V(q) = Vo, and of isotropic, parabolic va- 
lence and conduction bands that are perfectly nested in 
the absence of doping: 

k 2 

£kc = — £kr = Eq. (3) 

Ira 

Here Eq > 0, the case we consider, and Eq < imply 
semimetallic and semiconducting behavior, respectively. 
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Our results for ground states and excitations are 
based on unrestricted Hartree-Fock and time-dependent 
Hartree-Fock approximations respectively. Excitonic 
condensation is driven by interband exchange terms 
that favor spontaneous phase coherence. The intraband 
exchange self-energy has a spin-independent contribu- 
tion which we absorb into Eg and an important spin- 
dependent contribution which has not always been ac- 
counted for in previous work, but plays an essential role 
in stabilizing the ferromagnetic state JTo| . 

Dropping the direct Hartree interaction which plays no 
role for homogeneous states, the mean-field Hamiltonian 
is 
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Here 
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is the mean-field intraband exchange spin-splitting field 
(proportional to the band magnetization m Q ), 



Vr 
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are the singlet and triplet excitonic condensate order pa- 
rameters, and 



pZ><M) = 
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is the Hartree-Fock density matrix. 

We have solved these equations for a range of chemical 
potentials and external fields, allowing for density ma- 
trices that have spontaneous coherence between conduc- 
tion and valence bands and either spontaneous or induced 
spin-polarization in conduction and valence bands. Our 
results are summarized in Fig. [j]. There are in general 
several families of solutions, a normal state solution (N) 
in which no symmetries are broken, a ferromagnetic exci- 
tonic solution in which the excitonic condensate involves 
only conduction band electrons and valence band holes 
with particular spin orientations (see the inset in Fig. 
la)), and the excitonic insulator (EI) state in which both 
spin orientations participate in the excitonic condensate. 
Invariance of the Hamiltonian under independent spin- 
rotations of conduction and valence band systems in the 
absence of an external magnetic field is key to under- 
standing the phase diagram. For zero chemical potential 
and magnetic field a family M of symmetry-related EI 



states is lowest in energy. For h = 0, a family of ferro- 
magnetic states forms the ground state over a finite range 
of /i ^ 0. This family includes states with different spin 
orientations for condensed conduction and valence band 
particles, and states with different total magnetization 
magnitude. A weak external magnetic field selects from 
this family the state which has the largest magnetization 
magnitude aligned with the field direction. For /i > this 
is always a state in which spin down valence holes are con- 
densed. For large /z, where the largest contribution to the 
conduction band magnetization comes from uncondensed 
electrons, it follows that the collinear (COL) state which 
has coherence between parallel spin conduction and va- 
lence band electrons, is the ground state. For small , on 
the other hand, the non-collinear state (NC) is selected 
by a field. Both COL and NC ferromagnetic excitonic 
states compete with the EI state in which all spin orien- 
tations are paired, and with the N state, with the four 
different states separated by first order phase transitions. 
Both magnetic and pairing order parameters depend on 
both magnetic field and chemical potential, varying con- 
tinuously within a phase and discontinuously at a phase 
boundary. In region EI2 only spin down valence band 
electrons and spin up conduction band holes are con- 
densed. Unlike the NC state, this phase is insulating. 
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FIG. 1. Phase diagram in the presence of a field as a 
function of chemical potential [i a) or doping x b) and field 
h for Eg = 0.1 and Vo = 1.1. Total electronic density is 
given by no(l + x). Energy is measured in units of the band 
width 7r 2 /(2m) = 1. The adjacent pure phases occur as phase 
separated domains in the regions labelled D in b). The inset 
in a) shows a sketch of the Hartree-Fock quasiparticle bands in 
the ferromagnetic state. The inset in b) shows the dependence 
of doping on chemical potential. 



The slopes of the first order phase boundary lines are 
related to density and magnetization discontinuities by 
the following Clapeyron equation 



dh 
dfj, 
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Here x is the dimensionless quantity used to measure 
density and no is the density in the absence of dop- 
ing. As shown in the inset of Fig. lb) x always increases 
upon crossing these phase boundaries in the direction of 
increasing chemical potential. Positive slopes for these 
phase boundary lines correspond therefore to a magne- 
tization that decreases when crossing the phase bound- 
ary in the same direction. Inside each phase, density 
increases monotonically with chemical potential. It fol- 
lows that as a function of density, the phase diagram 
consists of pure phase regions interspersed with regions 
of phase separation. The regime consisting of phase sepa- 
rated COL and NC states, which have different magneti- 
zations and excitonic pairing of quite different character, 
is particulary unusual. Note also the reentrant normal 
state, predicted earlier O], induced in the COL by field 
at large x. 

Given the unusual nature of these ferromagnets, we 
can expect that their collective mode structure is un- 
usual. We have evaluated the dynamic linear response in 
a time-dependent Hartree-Fock approximation pOfl , We 
discuss here only spin-wave collective excitations, those 
that involve tilts of the magnetization orientation and 
are therefore off-diagonal in spin-indices when the valence 
and conduction band mean-field magnetization directions 
are chosen as quantization axes. 
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FIG. 2. Spin wave energies in units of the h = 0, fi = 
excitonic insulator order parameter Aq = 0.1 for the COL 
pure state. We measure the momentum in units of the lat- 
tice parameter. These results were obtained for E g — 0.1 
and Vq = 1.1, in units of the bandwidth. The main figures 
show zero field results, for two different values of the chem- 
ical potential corresponding to positive and zero conduction 
band magnetizations respectively The inset in b) shows the 
dependence of the COL state valence and conduction band 
magnetizations on fj,. Magnetization is measured in units of 
no. The inset in a) shows the collective excitation spectrum 
at finite field. 

Our results for the spin- wave collective modes E(q) 
of pure phase excitonic ferromagnets are summarized in 



Fig. |2|. For h — 0, there are two gapless collective modes 
with quadratic dispersion, corresponding to the indepen- 
dent conduction and valence band broken spin rotational 
symmetries, and an additional soft but gapped collec- 
tive mode. An interesting feature occurs in the doping 
dependence of this excitation spectrum. At an isolated 
value of \x the spectrum consists of one quadratically dis- 
persing mode and two linearly dispersing modes. To un- 
derstand the origin of this odd behavior, it is necessary 
to understand how the separate conduction and valence 
band contributions to the magnetization vary with dop- 
ing. Consider, for example, the collinear state. For the 
electron doping case illustrated here there are no un- 
condensed valence band holes. Uncondensed conduction 
band electrons, on the other hand, make a contribution 
to the magnetization opposite to that of the condensed 
conduction band electrons, yielding a net magnetization. 
At some point, the overall magnetization from the con- 
duction band vanishes. This is point at which the col- 
lective excitation spectrum is anomalous, and inciden- 
tally also the same point, at which COL and NC states 
have the same magnetization. Apparently the conduc- 
tion band excitation spectrum at this point is much like 
that of an antiferromagnet, with condensed and uncon- 
densed band electrons producing opposing magnetization 
contributions. 

In the presence of an external magnetic field our calcu- 
lations show that there are two Larmour collective modes 
at energy E = 2h = g^sH. At first sight the fact that 
conduction and valence collective modes both occur at 
this energy appears to suggest that their moments inter- 
act only with the external magnetic field and not with 
each other. If this were the case, however, each con- 
tribution to the magnetization would simply align with 
the field and there would be no energy barrier between 
collinear and non-collinear states. Furthermore, it is clear 
from our microscopic calculations that conduction and 
valence bands are coupled by the combined influences 
of excitonic pairing and Fermi statistics. This surpris- 
ing in-field collective mode behavior can be understood 
by considering the dependence of energy on the orien- 
tations of the conduction and valence band spins that 
participate in the pairing. In Fig. [| we plot, for both 
collinear and non-collinear states, energies obtained from 
self-consistent solutions of the mean-field equations for 
pairing between down spin valence band holes and con- 
duction band spins that have polar orientation angle 9 C 
with respect to the magnetic field. We see in Fig. || that 
there is indeed a non-trivial dependence of the energy 
on 9 C with the global minimum occurring at 9 C = n in 
the COL region of the phase diagram and at 9 C = 
in the NC region of the phase diagram. Why then do 
the conduction and valence band collective magnetiza- 
tion fluctuations not reflect this dependence on their rel- 
ative orientations? The answer again can be found in 
the symmetry of Hamiltonian. In a magnetic field, the 
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energy of the system is still invariant under independent 
rotations of conduction and valence band magnetizations 
around the field axis. Spin- wave collective excitations are 
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energy when magnetization is tilted away from the field 
direction. Since we are allowed to include only contribu- 
tions that are invariant under separate conduction and 
valence band spin-rotations, the quadratic coupling term 
rhx, c ' fn±.v, is not allowed. Higher order terms in the 
expansion of the magnetization, for example ones that 
go like (m± <c -rh± t c)('rh±,v ■fn±,v), and cause interactions 
between spin- waves, are symmetry allowed. 
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FIG. 3. Dependence of the energy on the angle between 
the spin of the condensed conduction band and the magnetic 
field for \i — 0.067 and fj, = 0.084 in which the ground state 
is at finite field in the NC and COL states respectively, and 
9 V = ty. From bottom to top h = 0.0, 0.01, 0.02. 

It is interesting to speculate on differences between 
the role of fluctuations in EI ferromagnets, compared to 
their role in other more familiar magnetic states. At 
low-temperatures, the Bloch T 3 / 2 M(T) law should still 
apply, despite the existence of several low-energy spin- 
wave branches. At higher temperatures, excitonic and 
magnetic order parameters vanish simultaneously at the 
mean-field level. We expect that mean-field theory criti- 
cal temperature estimates should be reliable in the weak- 
coupling limit Vo <C Eg. 

In closing, we remark that in any real excitonic 
insulator ferromagnet, corrections to the dominant 
term approximation would partially alter the picture 
outlined here. The largest corrections would likely 
come from interband interaction terms of the form 
^V cl. +tlC(T cl.,_ tiV(r/ Ci SL ' c(T 'C kv<T . This term will produces 
corrections to the Hartree-Fock energy of the form 
E HF oc 2|A S | 2 — m c m^ We find that the first term is 
always larger in magnitude, thus favoring triplet pair- 
ing and NC orientation in the ferromagnetic phase. The 
second term favors parallel conduction and valence band 
magnetizations at h = 0, coupling the two Goldstone 
modes and gapping one of them. The resulting h — /i 
phase diagram will depend on the competition of mag- 
netic energy and correction terms. In the weak coupling 
limit Ao <C \E g \, these terms are expected to be negligi- 
ble and the phase diagram qualitatively unchanged. 



[1] L.V. Keldysh and Yu. V. Kopaev, Sov. Phys. Solid State 

6, 2219 (1965). J. des Cloizeaux, J. Phys. Chem. Solids 

26, 259 (1965). 
[2] B.I. Halperin and T.M. Rice, Solid State Physics, 21, 

F. Seitz, D. Turnbull and H. Ehrenreich, eds., Academic 

Press, New York, 1968. 
[3] S. Conti, G. Vignale and A.H. MacDonald, Phys. Rev. B 

57, R6846 (1997) 
[4] LB. Spielman, J. P. Eisenstein, L.N. Pfeiffer, and K.W. 

West, Phys. Rev. Lett. 84, 5808 (2000);ibid 87,36803 

(2001). 

[5] K. Moon, H. Mori, K. Yand, S.M. Girvin, A.H. MacDon- 
ald, L. Zheng, D. Yoshioka and S.-C. Zhang, Phys. Rev. 
B 51, 5138 (1995). 

[6] B.A. Volkov, Yu. V. Kopaev and A.I. Rusinov., Sov. 
Phys. JETP 41, 952 (1975). B.A. Volkov, A.I. Rusinov 
and R. Kh. Timerov., Sov. Phys. JETP 43, 589 (1976) 

[7] D.P. Young, D. Hall, M.E. Torelli, Z. Fisk, J.L. Sarao, 
J.D. Thompson, H.-R. Ott, S. B. Oseeroff, R.G. Goodrich 
and R. Zysler, Nature 397, 412 (1999). 

[8] A. Hasegawa and A. Yanase, J. Phys. C 12, 5431 (1979). 
S. Massidda, A. Continenza, T.M. de Pascale and R. 
Monnier, Z. Phys. B 102, 83 (1997). 

[9] M.E. Zhitomirsky, T.M. Rice and V.I. Anisimov, Nature 
402, 251 (1999). 
[10] L. Balents and CM. Varma, Phys. Rev. Lett. 84, 1264 
(2000). MY. Veillette and L. Balents, Phys. Rev. B 65, 
14428 (2002). 

[11] V. Barzykin and L.P. Gor'kov, Phys. Rev. Lett. 84, 2207 
(2000). 

[12] M.E. Zhitomirsky and T.M. Rice, Phys. Rev. B 62, 1492 
(2000). 

[13] L. Balents, Phys. Rev. B 62, 2346 (2000). 
[14] T. Ichinomiya, Phys. Rev. B 63, 5113 (2001) 

M 



[16 
[17 



[18 
[19 

[20 



0010344 



T.A. Glo or, M.E. Zhitomirsky and T.M. Rice, |cond-mat 
R. Shindou, N. Nagaosa and A.S. 



S. Muraka mi 

Mishchenko, |cond-mat 0107215 . 
J.D. Denlinger, J.A. Clack, J.W. Allen, 
D.M. Poirier, C.G. Olson, J.L. Sarrao, Z. Fisk, |cond-mat 



G.-H. Gweon, 



0009022, i 



: 0107429 
P. van Gelderen, 



H.J. Tromp, P.van Gelderen, P.J. Kelly, G. Brocks and 
P.A. Bobbert, Phys. Rev. Lett 87, 16401 (2001) 
Some past work has incorrectly predicted energy differ- 
ences between states that are related by transformations 
that are in the symmetry group of the Hamiltonian. 
In the EI case we recover the acoustic modes previ- 
ously predicted in A.N. Kozlov and L.A. Maksimov, Sov. 
Phys. JETP 22, 889 (1966) and D. Jerome, T. M. Rice 
and W. Kohn, Phys. Rev. 158, 462 (1967). Due to the 



4 



SU(2) x SU(2) symmetry of the dominant term approxi- 
mation, this branch is 4-degenerate. Two of these modes 
corresponds to the phonons of the condensates while the 
character of the other two depends on the properties of 
the selected ground state, e.g. the magnons of the spin 
density wave in the triplet state. 
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